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Intestinal immunoglobulin A (IgA) ensures host defense and symbiosis with our commensal 
microbiota. Yet previous studies hint at a surprisingly low diversity of intestinal IgA, and it is 
unknown to what extent the diverse Ig arsenal generated by somatic recombination and diversi- 
fication is actually used. In this study, we analyze more than one million mouse IgA sequences 
to describe the shaping of the intestinal IgA repertoire, its determinants, and stability over time. 
We show that expanded and infrequent clones combine to form highly diverse polyclonal IgA 
repertoires with very little overlap between individual mice. Selective homing allows expanded 
clones to evenly seed the small but not large intestine. Repertoire diversity increases during 
aging in a dual process. On the one hand, microbiota-, T cell-, and transcription factor ROR-yt- 
dependent but Peyer's patch-independent somatic mutations drive the diversification of 
expanded clones, and on the other hand, new clones are introduced into the repertoire of aged 
mice. An individual's IgA repertoire is stable and recalled after plasma cell depletion, which is 
indicative of functional memory. These data provide a conceptual framework to understand the 
dynamic changes in the IgA repertoires to match environmental and intrinsic stimuli. 



The intestine represents the major Ig-producing 
compartment in mammals (Brandtzaeg and 
Johansen, 2005; Macpherson et al., 2008). 
Most intestinal plasma cells produce IgA and 
express the J chain, allowing for transepithelial 
transport and release of secretory IgA (SIgA) 
into the gut lumen (Brandtzaeg, 2009). Among 
the many functions of SIgA are the vaccination- 
induced protection against cholera toxin (Lycke 
and Holmgren, 1986) and protection against 
Salmonella enterica (Wijburg et al., 2006) and 
influenza infection (Asahi et al., 2002). In addi- 
tion to these classical immunoprotective func- 
tions, IgA plays a central role in regulating 
host-microbiota homeostasis. Few plasma cells 
are present in neonates and in germ-free mice, 
but plasma cell numbers rapidly increase upon 
colonization with commensal microbiota, and 
even transient exposure to live bacteria induces 
a stable plasma cell population (Hapfelmeier 
et al., 2010). Colonization of germ-free mice 
lacking B cells results in enhanced translocation 
of bacteria to gut-draining mesenteric LNs 
(mLNs) compared with WT mice, and SIgA 
deficiency increases titers of flora-binding IgG 
in serum (Johansen et al., 1999; Macpherson 
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and Uhr, 2004; Sait et al., 2007). These obser- 
vations show that SIgA limits the access of in- 
testinal microbiota to systemic tissues. However, 
mechanistically, this process is poorly understood, 
and even less is known about the specificity 
and diversity of the IgA repertoire. Intestinal IgA 
binds to the commensal microflora, and such 
reactivity seems to be specific for distinct bacterial 
epitopes (Hapfelmeier et al., 2010). Moreover, 
screening a set of IgA hybridoma lines, no more 
than ^20% of antibodies were rated polyreactive, 
whereas most clones did not show cross-reactivity 
(Benckert et al., 2011). This hints at a diverse 
IgA repertoire, potentially reflecting the com- 
plexity of the intestinal microbiota and anti- 
genic load. Once induced, microbe-specific IgA 
levels are comparably stable over time but be- 
come attenuated when new bacterial species are 
introduced to the microbiota (Hapfelmeier 
et al., 2010). Therefore, the IgA repertoire might 
be dynamically shaped to mirror the actual 
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composition of the intestinal microbiota. Because changes in 
microbiotal composition are rapid during infancy and early life 
(Dominguez-Bello et aL, 2011), one might speculate that the 
IgA repertoire could undergo similarly profound remodeHng 
during this period of life. Moreover, mice capable of SIgA pro- 
duction and secretion but with an impaired capacity to acquire 
somatic mutations (SMs) showed a dysbiosis of their microbiota 
and hyperplasia of gut-associated lymphoid tissues (GALTs; Wei 
et al., 2011). This indicates that not only do microbiota induce 
specific IgA responses, but the specificity of the IgA repertoire 
may in turn also affect the composition of the microbiota. 

Considering the near-infinite number of potential CDR3 
(complementarity determining region 3) sequences and the 
complexity of the antigen load, it is surprising that the study of 
IgA repertoire diversity by CDR3 length analysis and se- 
quencing hints at an oligoclonal IgA repertoire of low diver- 
sity. In particular, the frequent observation of clonally related 
plasma cells (Dunn-Walters et al, 1997; Holtmeier et al, 2000; 
Stoel et al, 2005, 2008; Yuvaraj et al, 2009) was interpreted 
to indicate low repertoire diversity. However, previous studies 
used conventional sequencing technologies and had to ex- 
trapolate from the analysis of few sequences to the overall IgA 
repertoire. In this study, we used high-throughput sequencing 
to investigate the shaping of the IgA repertoire, its determi- 
nants, and stability over time, as well as differences between 
individuals. Analyzing more than one million V H sequences, 
we show that the IgA repertoire comprised both highly 
expanded and low frequency clones. These two components 
combined to make a highly polyclonal IgA repertoire. Ex- 
panded clones were evenly distributed along the small but not 
large intestine and dominated the repertoire in young mice. 
However, during aging, repertoire diversity increased by the 
ongoing accumulation of low frequent clones and microbiota-, 
T cell—, and transcription factor ROR7t— dependent but 
Peyer's patch (PP)— independent hypermutation. Moreover, 
expanded clones were rapidly recalled after plasma cell de- 
pletion, indicating that the intestinal IgA system possesses 
memory characteristics. 

RESULTS 

Expanded and low frequency clones constitute 
polyclonal Individual IgA repertoires 

To establish comprehensive IgA repertoire analysis, RNA was 
isolated from whole small intestine (SI) tissue and transcribed 
into cDNA, and V H sequences were amplified by PCR with 
primers annealing to FR1 (frame work region 1) and the Ca 
domain. Amplicons of ~450-bp length were sequenced by 
454 sequencing, and for each sample 5,000—30,000 sequences 
of expected length, including the flanking Ca domain primer 
and encoding intact Ig, were analyzed. Previous studies described 
clonally related sequences even when comparably few sequences 
were analyzed (Dunn- Walters et al., 1997, 2000; Holtmeier 
et al., 2000; Stoel et al, 2005; Yuvaraj et al, 2009). Indeed, in 
99.7 + 0.3%, we observed at least one identical sequence pair 
within a set of 25 CDR3 sequences (5,000 randomly picked 
sets of 25 sequences obtained from six 10— 11-wk-old mice), 



and some clones accounted for >1% of all sequences. However, 
along with such highly expanded clones, we observed unique 
sequences. In a set of 5,000 sequences analyzed per mouse, 
25-35% of all full-length V H sequences and 3.9 ± 0.8% of all 
CDR3 sequences were obtained only once (mean + SD in six 
10— 11-wk-old C57BL/6 mice). This suggests that the IgA 
plasma cell pool includes highly expanded clones along with low 
frequency clones. Therefore, the presence of clonally related 
sequences observed in small sequence sets does not by itself indi- 
cate low diversity and oligoclonality of the entire repertoire. 

To visualize this feature of the IgA repertoire, we sorted all 
different CDR3 sequences in a set of 5,000 sequences from most 
frequent to single sequences. These sorted sequence sets were 
aligned along the x axis, and on the y axis we indicated how 
many times the respective sequence was observed (Fig. 1 A). This 
diagram illustrates the distribution of the IgA repertoire, which 
is composed of expanded and lowly frequent clones. To allow 
for further analysis of the data, each CDR3 sequence was 
assigned a rank based on its abundance. Sequences obtained 
only once were given rank 1, sequences obtained twice 
rank 2, and high abundant sequences given a rank equaling 
to how many times the respective sequence was observed. We 
next detemiined the frequency of each rank, e.g., among a rep- 
resentative set of 5,000 CDR3 sequences, we identified 294 
sequences of rank 1; thus, the frequency of rank 1 was set to 294. 
Log transformation of rank and frequency allowed formally 
separating two components of the repertoire. A first compo- 
nent followed the power law, i.e., appeared as straight line, 
whereas the second component presented as a long tail in log 
frequency-rank diagrams (Fig. 1 B). In this respect, the IgA 
repertoire resembled previously described distributions of 
human CD8 + memory T cells (Naumov et al., 2003, 2011) and 
mouse regulatory T cells (Haribhai et al., 2011). Calculating 
the cutoff point separating both components (Naumov et al., 
2003) showed that in a set of 5,000 sequences, sequences ob- 
tained less than roughly 20 times (19.4 + 3.2; n = 6) contrib- 
uted to the first component, whereas sequences obtained more 
frequently belonged to the second component of the repertoire. 
Notably, expanded clones contributed substantially to the 
total sequence pool but made up only a small proportion of all 
different CDR3 sequences (Fig. 1 C). 

The repertoire's two components might constitute an in- 
trinsic feature of the plasma cell repertoire. Alternatively, the 
composite nature of the IgA repertoire observed in our whole 
tissue samples might reflect a mixture of different cell types 
carrying largely divergent amounts of IgA encoding messenger 
RNA, such as bona fide plasma cells expressing high levels of 
IgA transcripts and nonplasma cells expressing much lower 
amounts of IgA encoding messenger RNA. To distinguish be- 
tween these alternatives, we analyzed IgA + CD138 + plasma 
cells sorted from the SI lamina propria (Fig. 1 D). Analyzing 
the sequence repertoire of sorted IgA + CD138 + cells, we ob- 
served, similar to our observation obtained from whole tissue, 
highly expanded along with low frequent sequences (Fig. 1 E). 
Moreover, log transformation of rank and frequency revealed 
two distinct components formally separated by a cutoff point 
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similar to what we observed analyzing whole tissue samples 
(cutoff point for sorted plasma cell populations, 23.9 + 3.4; 
n = 4; Fig. 1 E). We conclude that the two components are a 
fundamental characteristic of the intestinal IgA repertoire that 
needs to be considered when the diversity and nature of the 
IgA repertoire is discussed. Previous work used measures of 
ecosystem diversity, which considers the richness and evenness 
of a population, to describe the diversity of T cell receptor 
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Figure 1. The IgA repertoire comprises highly expanded and low 
frequency clones. (A) In a set of 5,000 sequences, all different CDR3 se- 
quences were listed along the x axis from frequent to single sequences, and 
the number of sequence reads for each different CDR3 sequence were dis- 
played on the y axis. Each slice represents an independently analyzed 
mouse. (B) Log(frequency) versus log(rank) diagram for one representative 
mouse reveals two components of the CDR3 population. A first component 
(I) follows power law characteristics and comprises clones present at low 
frequency. A second component (II) comprises highly expanded clones and 
appears as a long tail in the diagram and is marked in red. (C) Bars depict 
the contribution of both components to all CDR3 sequences and different 
CDR3 sequences (mean + SD). All data are based on the analysis of six 10- 
1 1-wk-old C57BL/6 mice housed under SPF conditions. (D) Representative 
pseudocolor plot demonstrating IgA/CD 1 38 staining gated on live (DAPh) 
single SI lamina propria cells. (E) Log (frequency) versus log(rank) diagram of 
5,000 CDR3 sequences obtained from sorted lgA + CD138 + plasma cells 
pooled from two 10-wk-old C57BL/6 mice. Results are representative of 
four independent experiments performed. 



repertoires (Naumova et al., 2009). In the case of our data, the 
richness of the IgA population correlates to the number of 
different CDR3 sequences and positively correlates with di- 
versity. Because lowly frequent clones account for the majority 
of all different CDR3 sequences, the repertoire's nonexpanded 
first component is particularly relevant when the richness of 
the IgA repertoire is estimated. We therefore conclude that in 
contrast to previous suggestions (Stoel et al., 2005; Yuvaraj 
et al., 2009), the intestinal IgA repertoire is highly polyclonal. 

We next compared the IgA repertoire between three pairs 
of C57BL/ 6 mice housed under specified pathogen— free (SPF) 
conditions in three different cages. In all mice, consistently the 
majority of V H domain sequences belonged to the V H 1 or V H 5 
family (48.1 ± 13.8% V H 1 and 32.1 ± 14.3% V H 5; n = 6). 
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Figure 2. Individual mice show largely nonoverlapping IgA reper- 
toires in the intestine. (A) In a set of 5,000 sequences per sample, all 
VDJ combinations were enumerated, and the sequence sets were sorted 
from most frequent to least frequent as observed in mouse number 1, 3, 
or 5 (indicated in red). The number of sequence reads assigned to each 
VDJ combination is indicated by color code. Similarly, sequence sets were 
sorted according to the frequencies of CDR3 amino acid sequences. 
Please note that logarithmic transformation was used to generate the 
color code and even sequences observed only three times are nonblack. 
(B) CDR3 amino acid and nucleotide sequences were compared between 
different mice, and the number of identical sequences was enumerated. 
The number of different CDR3 sequences in a given animal appears in the 
diagonal of the table. Very few CDR3 sequences were shared between 
different animals. All data are based on the analysis of six 10-1 1-wk-old 
C57BL/6 mice housed under SPF conditions. 
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Figure 3. Identical CDR3 sequences seed different fragments of 

the SI. (A) Nucleotide CDR3 sequences were compared between proximal 
jejunum (SM), distal jejunum (SI2), ileum (SI3), and colon (c) of 13-wk-old 
C57BL/6 mice, and the number of identical sequences was enumerated. 



Analyzing the repertoire's second component, i.e., highly 
expanded sequences, revealed that some VDJ combinations 
were abundantly used throughout several individual mice 
(Fig. 2 A). Still, when comparing CDR3 amino acid sequences 
in different mice, only very few CDR3 sequences expanded 
in one mouse were also detected in another animal (Fig. 2, 
A and B). Specifically, only three CDR3 sequence shared 
between two mice were highly expanded in both individuals, 
and no CDR3 amino acid sequence was observed in all mice. 
We quantified this observation by calculating Morisita-Horn 
indices (MHIs; Magurran, 1988), which score two identical 
populations as 1 and two completely disparate populations as 0. 
Comparison of CDR3 sequence repertoires of different mice 
consistendy resulted in MHIs close to 0 (MHI 0.0013 ± 
0.0020; n = 6), which is indicative of dissimilar populations. 
This indicates that even when exposed to the identical envi- 
ronment and possibly sharing a similar composition of micro- 
biota, there is little conversion of the CDR3 sequence 
repertoire between different mice. 

CCR9 differentially targets expanded clones to the small 
and large intestine 

Although the sequence repertoire of different mice showed 
litde overlap, we readily detected identical CDR3 sequences 
in different parts of the intestine obtained from the same ani- 
mal (Fig. 3 A). CDR3 sequences from different SI fragments 
from the same mouse were significantly more similar than 
SI and colon (Fig. 3 B). Still, the sequence repertoire in SI and 
colon of the same mouse was more similar compared with 
repertoires obtained from different mice (Fig. 3 B). Notably, 
largely overlapping sequence repertoires between indepen- 
dently processed samples of the same but not different mice 
also ensured that expanded sequences did not result from 
a PCR bias and that sequence errors did not preclude conclu- 
sive data analysis. Overlap of different SI samples resulted 



Data are based on 10,000 sequences per sample, and the number of dif- 
ferent CDR3 sequences in each sample appears in the diagonal of the 
table. (B) MHIs were calculated comparing CDR3 sequence pools obtained 
from various gut fragments in WT and CCR9~'~ mice (WT, n = 2; 
CCR9~'~, n = 2). Individual mice are indicated by open and closed sym- 
bols, and comparisons calculated are SI1 to SI2, SI1 to SI3, and SI2 to SI3 
to describe similarity of SI fragments and SI 1 to c, SI2 to c, and SI3 to c 
for comparison of SI- and large intestine-derived sequence pools. Com- 
parison of CDR3 sequence pools between different mice yields MHIs close 
to 0 (first and second columns). Horizontal lines indicate the mean. 
*, P < 0.05;**, P < 0.001. (C) Venn diagrams illustrate the overlap of CDR3 
sequence pools from SI1, SI3, and colon in one WTand one CCR9~'~ mouse. 
Numbers indicate overlap in percentages. (D) 10,000 CDR3 sequences 
from SM, SI2, SI3, and colon were sorted according to their frequency in 
SI1 (WT, left) or colon (WT, right). Because overlap of sequence sets is 
driven by expanded sequences, we limited our analysis to the 100 most 
frequent sequences and visualized their frequency (in percentages of all 
sequences) by color codes. Note that upon sorting according to frequency 
in the proximal jejunum, more sequences with similarly high frequency 
could be found in CCR9~'~ mice than in WT mice. 
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primarily from overlapping expanded clones (Fig. 3 C), i.e., 
the repertoire's second component. Sorting 10,000 CDR3 
sequences obtained from the proximal jejunum from most 
frequent sequences to sequences observed only once, we found 
that the identical CDR3 sequences were similarly frequent 
in distal jejunum and ileum (Fig. 3D). In contrast, many se- 
quences abundantly present in the SI were infrequent or unde- 
tectable in the colon, and only a minority of sequences was 
present at high frequency in both small and large intestinal 
sequence sets. Sorting the sequence sets, according to their fre- 
quency in the colon, we found a comparable scenario, i.e., 
only a minority of the expanded sequences dominating the 
sequence repertoire in the colon was similarly frequent in 
the SI (Fig. 3 D). 

Differences in the plasma cell repertoire between SI and 
colon might be caused by selective homing of plasmablasts or 
alternatively stem from subsequent local events in the lamina 
propria. To study the role of plasmablast homing, we used 
mice lacking CCR9 (chemokine [CC motif] receptor 9), 
which mediates homing of plasmablasts into the SI but not 
colon (Pabst et al., 2004). In CCR9 _/ ~ mice, SI- and large 
intestine— derived sequence pools showed higher similarity 
(Fig. 3 B) and larger sequence overlap (Fig. 3, C and D) as 
compared with WT mice. Because CCR9 is down-regulated 
on plasma cells inside the lamina propria (Pabst et al., 2004), 
these data indicate that expansion of distinct specificities 
might occur before plasma cell homing to the lamina propria. 
Yet CCR9 expression is not restricted to plasmablasts, and 
we cannot rule out the possibility that other CCR9-expressing 
cells contribute to the increased repertoire overlap observed 
in CCR9~ 7 ~ mice. To directly examine the extent of plasma 
cell proliferation inside the lamina propria, we compared the 
distribution of recently proliferated plasma cells after 7 d of 
continuous BrdU application with distributions predicted in 
silico. In a representative experiment, we observed a total of 
180 BrdU + IgA + cells in 134 villi (pooled from five mice). 
According to this experimental observation, we predicted 
distributions of 180 cells throughout 134 villi in silico assum- 
ing that cells did not proliferate, performed one division, or 
performed two divisions (Fig. 4 A). In silico simulations based 
on two divisions were distinguished by the noticeable accu- 
mulation of numerous cells within few villi, whereas 69% of 
all villi were devoid of any cell. Similarly, distributions ob- 
tained for single divisions inside the gut predicted that ^49% 
of all villi would not contain any cell (Fig. 4 A). Both of these 
predicted properties did not resemble the experimentally 
observed spreading of BrdU + IgA + cells throughout individual 
villi. In contrast, the experimental data, i.e., the observed num- 
ber of BrdU + IgA + per villus, fitted simulations that did not 
assume in situ proliferation (Fig. 4 B). Yet a similarly good 
representation of the experimental data, i.e., similarly low sum of 
squared errors, was obtained for models that assumed that 10% 
or less of all cells divided once inside the lamina propria. We 
therefore conclude that plasmablasts evenly homed to the SI but 
not colon, and in situ proliferation seems unlikely to fully ex- 
plain the presence of expanded clones in the IgA repertoire. 



Age, microbiota, and T cells drive repertoire diversity 

Besides CDR3 diversity, somatic hypermutation (SHM) 
contributes to IgA diversity, and in humans most IgA se- 
quences are highly mutated (Barone et al., 2011; Gibbons 
and Spencer, 2011). Thus, we extended our analysis from 
CDR3 to full-length V H sequences and frequencies of SMs 
to investigate further determinants of repertoire diversity. In 
a first set of experiments, we compared the IgA repertoire 
in mice of various age. We did not observe significant dif- 
ferences with respect to V H family usage between mice of 
various age, but the frequency of CDR3 sequences lacking 
N nucleotides (N~ sequences) was higher in 4-wk-old mice 
compared with older mice (44.6 + 6.2 in 4-wk-old [n = 2] 
compared with 11.2 + 9 in 10— 11-wk-old mice [n = 6]). 
Lack of nontemplated N nucleotides is a characteristic feature 
of peritoneal B cells and B cells of fetal origin (Vale et al., 
2010). Thus, we speculate that a higher frequency of N~ 
sequences might result from B cells generated in the fetal 
liver prevailing in the repertoire of young mice. Yet these 
differences did not results in more similar IgA repertoires in 
young mice compared with old mice, and also comparison of 




1 134 



Figure 4. Newly generated plasma cells are evenly distributed in 
the gut mucosa. (A) A random distribution of 180 cells throughout 134 
individual villi was predicted by in silico simulation assuming no (black), 
one (green), or two (blue) cell divisions. (B) Mice received BrdU continu- 
ously with the drinking water for 7 d, and the number of BrdU + lgA + cells 
in individual villi was enumerated by immunofluorescence microscopy. 
Inspecting 134 villi (pooled from five mice), we found 1 80 BrdU + lgA + 
cells. The x axis corresponds to individual villi sorted from villi containing 
no cells to villi containing multiple BrdU + lgA + cells, and each villus is 
represented by a vertical red line. Mean and 95% confidence interval 
of the in silico simulation with no cell divisions are depicted by a black 
line and blue-filled area. Similar results were observed after 2 d 
of BrdU application. 
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CDR3 sequence repertoires of 4— 5-wk-old mice resulted in 
an MHI close to 0 (MHI 0.00442 ± 0.00461; n = 4). Differendy 
aged mice showed striking differences in their repertoires 
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Figure 5. IgA repertoire diversity increases with age. (A) Diagrams 
depict the number of SMs observed in representative 4-, 6-, 10-, and 18- 
wk-old mice. (B) Bars depict mean + SD SM frequencies in FR1, FR2, FR3, 
CDR1, and CDR2 observed in the following number of animals: 4-5 wk, 
n = 4; 6 wk, n = 2; 10-1 1 wk, n = 6; and 18-22 wk, n = 5. (C) CDR3 
sequences were sorted by rank, from high to low frequency, and the number 
of identical sequences was plotted (compare with Fig. 1 A). Please note 
that the number of unique sequences increased with age. (D) Bars depict 
the contribution of both components to all CDR3 sequences (mean + SD). 



when we analyzed frequencies of SMs and CDR.3 distribu- 
tions. Whereas 4— 5-wk-old mice showed few SMs, their fre- 
quencies gradually increased in 6-, 10—11-, and 18— 22-wk-old 
mice (Fig. 5, A and B). In 18— 22-wk-old WT mice housed 
under SPF conditions, most V H sequences showed numerous 
SMs. SM frequencies were highest in CDR1 and CDR2 
regions, and the number of replacement mutations exceeded 
that of silent mutations (Fig. 5 B) . This indicates that age is 
a critical determinant of SM frequencies. Moreover, age also 
affected CDR3 sequence distributions. In particular, ex- 
panded clones, i.e., the second component not following 
power law characteristics, were more abundant in young 
compared with old mice, and conversely, aged mice accumu- 
lated increasing numbers of CDR3 sequences present at low 
frequency (Fig. 5 C). The increase in lowly frequent clones 
in aged mice is also reflected in a shift of the repertoire's two 
components (Fig. 5 D). This shift may indicate the ongoing 
recruitment of low frequent sequences into the repertoire. 
Yet alternatively low SM frequencies and higher frequency 
of N~ sequences in young compared with aged mice may 
contribute to observed differences. Notably, expanded clones 
present in 18— 22-wk-old mice contained as many SMs as low 
frequency clones (Fig. 6 A), and expanded clones of identical 
CDR3 showed complex phylograms, which frequently con- 
tained several major branches (Fig. 6 B). We suggest that 
during aging, two independent processes contribute to the 
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Figure 6. Expanded clones accumulated extensive SMs. (A) Dia- 
grams depict the frequencies of SMs in two representative 18-wk-old 
mice (out of five 18-22-wk-old mice analyzed) in sequences comprising 
the 10 most frequent CDR3 sequences or in V H sequences carrying CDR3 
sequences that were present only once in a set of 17,000 (left) and 23,000 
(right) sequences. (B) Phylogram demonstrating clonal relationship of the 
most frequent CDR3 sequences in one representative 11- and 18-wk-old 
mouse. Numbers indicate the number of SMs in selected sequences. 
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Figure 7. SHM requires T cells, microbiota, and ROfVyt but not 
PPs. (A) Germ-free mice were colonized by cohousing with sentinels 
raised under SPF conditions for the indicated number of days. Symbols 
indicate the total number of lgA + and Ig M + cells in individual mice 



diversification of the repertoire. First, on expense of expanded 
clones dominating the repertoire in young mice, the number 
of CDR3 sequences observed at low frequency increased in 
older mice, and second, SMs accumulated. 

Next, we analyzed the role of microbiota in generating 
repertoire diversity. Whereas germ-free mice harbored only 
few intestinal plasma cells, cohousing of germ-free mice with 
SPF mice resulted in rapid induction of plasma cells in the 
lamina propria. After 4 wk, the number of IgA + and IgM + 
cells was indistinguishable between ex— germ-free mice and 
sentinels housed under SPF conditions (Fig. 7 A). Notably, 
IgA + and IgM + cell numbers increased with similar kinetics 
during the filling of the lamina propria in ex— germ-free mice 
(Fig. 7 A). This indicates that most intestinal IgA-secreting 
cells undergo class switch recombination before entering 
the lamina propria and is consistent with the lack of overt 
proliferating plasma cells in the lamina propria described in 
the previous section. The CDR3 sequence repertoire in 
germ-free mice showed the characteristic two components, 
low frequent along with highly expanded clones, also ob- 
served in SPF-reared mice (Fig. 7 B), indicating that a polyclonal 
IgA repertoire can be generated independent of microbial 
stimulation. Still, germ-free mice showed few SMs compared 
with age-matched SPF mice (Fig. 7, C and D). Even though 
after 4 wk of cohousing plasma cell numbers were indistin- 
guishable, V H domain sequences from ex— germ-free mice 
still showed fewer SMs than those from age-matched SPF 
mice (Fig. 7, C and D), indicating that SHM was delayed 
compared with accumulation of plasma cells. Similarly, in 
T cell— deficient mice, we observed a polyclonal CDR3 pop- 
ulation, whereas SM frequencies were rather low. Notably, 
mice treated in utero with lymphotoxin-(3-receptor fusion 
protein, which blocks the development of PPs, showed 
similarly high SM frequencies compared with PP-sufficient 
mice (Fig. 7, C and D). Besides PPs, the SI harbors a large 
number of small-sized lymphoid aggregates, which in response 
to microbial stimulation can develop into B cell— rich iso- 
lated lymphoid follicles (ILFs; Pabst et al., 2006). Because, 



observed by histology counting six fields of view from at least three dif- 
ferent sections per mouse covering different fragments of the SI. Hori- 
zontal lines indicate the mean. (B) In a set of 5,000 sequences, all 
different CDR3 sequences were listed along the x axis from frequent to 
single sequences, and the number of sequence reads for each different 
CDR3 sequence was displayed on the y axis. Each slice represents an inde- 
pendently performed experiment as indicated. (C) Diagrams depict the 
number of SMs observed in representative 10-wk-old WT (see Fig. 5), 
germ-free, ex-germ-free, CD3~'~, PP-free, and R0Fi7t~'~ mice. 
(D) Bars depict mean + SD SM frequencies in FR1, FR2, FR3, CDR1, and 
CDR2 observed in the following number of animals: 10-1 1-wk-old WT 
mice (n = 6; identical to data in Fig. 5 B and shown to improved 
readability), germ-free mice (9 wk; n = 2), ex-germ-free mice (1 1 wk; 
n = 2), CD3~'~ mice (10 wk; n = 3), mice lacking PPs (12 wk; n = 2), 
and R0R-/t~'~ mice (20 wk; n = 2). For IgA sequences from 
PP-free mice, RNA was isolated out of sorted lgA + CD138 + intestinal 
plasma cells. For each mouse, a minimum of 10,000 sequences 
were analyzed. 
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first, ILFs have been suggested to support PPs and T cell- 
independent IgA switch (Tsuji et al., 2008) and, second, our 
experimental approach of PP ablation results in increased 
formation of ILFs (McDonald et al., 2005), we next analyzed 
the IgA repertoire in ROR7t~ / ~ mice. ROR"yt _/_ mice 
lack lymphoid tissue inducer cells and in consequence do not 
develop LNs, PPs, and ILFs (Eberl and Littman, 2004). Still, 
aged ROR7t~ / ~ mice generate intestinal plasma cells (Tsuji 
et al., 2008). In contrast to PP-deficient mice, 20-wk-old 
ROR~yt~ /_ mice showed very few SMs, similar to the situa- 
tion in T cell— deficient mice (Fig. 7, C and D). No differ- 
ences were observed in any of these mice with respect to V H 
family usage. In conclusion, these findings indicate that 
CDR3 diversity could be generated independent of micro- 
biota, T cells, and GALTs. In contrast, diversification of the 
repertoire by SHM required microbiota and T cells but not 
PPs. Still, SMs were absent in aged ROR7t~ / ~ mice, even 
though these mice possess T cells. This may indicate that 
ROR~yt-dependent lymphoid tissue inducer cells might di- 
rectly contribute to SHM. Alternatively, lack of SMs in 
ROR~yt~ /_ mice might be a consequence of impaired lym- 
phoid organogenesis. We therefore propose that different 
lymphoid compartments, potentially comprising PPs, mLNs, 
ILFs, and spleen, might have overlapping functions in driving 
T cell— dependent accumulation of SMs. 

The IgA system recalls previously selected specificities 
after plasma cell depletion 

To investigate the stability of individual repertoires, we es- 
tablished an experimental setup that allows on the one hand 
to deplete intestinal plasma cells and on the other hand to 
follow an individual's repertoire over time. To deplete intes- 
tinal plasma cells, we used the proteasome inhibitor Bortezo- 
mib, which is used for the treatment of myeloma patients. 
Bortezomib has been shown to deplete plasma cells in bone 
marrow and spleen (Neubert et al., 2008) and has less severe 
effects on other immune cells. We observed that two consec- 
utive injections of Bortezomib depleted most intestinal IgA- 
secreting cells (Fig. 8, A and C). Yet numbers of IgA-secreting 
cells rapidly recovered after depletion, and after 5—7 d, plasma 
cell numbers in the SI had returned to normal values (Fig. 8, 
A and C). Virtually all of these plasma cells were newly gener- 
ated (Fig. 8, B and C). Replenishment of plasma cell numbers 
after depletion showed similar kinetics in PP-deficient and 
splenectomized mice compared with WT mice (Fig. 8 D). 
However, plasma cell repopulation was delayed in CCR9 _/_ 
mice (Fig. 8 D). This suggests that the plasma cell population 
installed after depletion was derived from circulating progen- 
itors that required CCR9 to enter the lamina propria. Fur- 
thermore, the plasma cell repertoire reinstalled after depletion 
showed SM frequencies similar to nondepleted age-matched 
controls (Fig. 9 A). This indicates that the plasma cell pool 
was reestablished from a previously selected and mutated cell 
population. Moreover, the reinstalled IgA repertoire showed 
the characteristic two components (Fig. 9 B) without major 
changes in their contribution (Fig. 9 C). To compare the 
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Figure 8. Intestinal IgA-secreting plasma cell numbers rapidly recover 
after depletion. (A) Mice were injected i.p. on two consecutive days with Bor- 
tezomib (Bz), lamina propria cells were isolated, and IgA-secreting plasma cells 
were quantified by ELISP0T on days 0, 2, 4, and 7 as indicated. Circles represent 
individual mice analyzed and pooled from two or more experiments, and hori- 
zontal lines indicate the mean. (B) Mice received BrdU with the drinking water 
starting 6 d before Bortezomib injection as indicated. Lamina propria cells were 
isolated, and the frequency of BrdU + cells among CD138 + lgA + cells was deter- 
mined by flow cytometry. Control mice did not receive BrdU. (C) SI tissue was 
embedded in paraffin, sectioned, and immunostained for IgA (green) and BrdU 
(red). Images are representative of at least five mice analyzed in two independent 
experiments. Boxed areas are shown at higher magnification in the bottom pan- 
els. Bars: (top) 50 pm; (bottom) 10 pm. (D) Plasma cell numbers rapidly recover 
after depletion in PP-deficient mice and splenectomized mice but not in CCR9~'~ 
mice. Mice received two consecutive injections of Bortezomib, and plasma cell 
numbers were quantified by ELISP0T in nondepleted mice and at days 0, 5, or 7 
after depletion as depicted. Symbols represent individual mice pooled from at 
least two experiments. Horizontal lines indicate the mean. *, P < 0.05. 
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Figure 9. Identical specificities were recalled after plasma cell 
depletion. (A) Bars depict mean + SD SM frequencies in FR1, FR2, FR3, 
CDR1, and CDR2 present in mice after plasma cell depletion or without 
depletion (four to five mice per group). (B) Representative log(frequency) 
versus log(rank) diagram of 5,000 CDR3 sequences. (C) Bars depict the 
contribution of both components to all CDR3 sequences (mean + SD). 
(D) 10,000 CDR3 sequences were compared between SI biopsies (biop) 
and SI of two control and two Bortezomib (Bz)-treated mice, and the 
number of identical CDR3 sequences was enumerated. Overlap of biopsy 
and SI from the same mouse is highlighted by red numbers. (E) MHIs were 
calculated comparing CDR3 sequences from the biopsies with the corre- 
sponding SI. Mice were 11-15 wk old at the time of sacrifice. Horizontal 
lines indicate the mean. (F) Venn diagrams illustrate percent overlap of 
the 100 most frequent CDR3 sequences and 900 randomly picked unique 
sequences from the SI, biopsy, and spleen of one untreated mouse. 



plasma cell repertoire before depletion and after plasma cell 
replenishment within the same animal, we next obtained SI 
biopsies before Bortezomib injection. Interestingly, the rep- 
ertoire observed in biopsies was similar to the repertoire 
recalled after plasma cell depletion (Fig. 9, D and E). More- 
over, we noted a substantial overlap of the IgA repertoire in 
the SI and spleen of the same animal (Fig. 9 F). This indicates 
that after plasma cell depletion, the IgA system recalled previ- 
ously selected specificities. In sum, our results suggest that an 
individual's IgA repertoire undergoes ongoing modulation, 
by recruitment of new nonexpanded clones into the popula- 
tion and on the basis of already established expanded clones. 
Diversification of expanded clones might require the repeated 
entry and selection of circulating precursors in germinal cen- 
ters, which might also underlie the striking stability of an 
individual's IgA repertoire. 

DISCUSSION 

The intestinal mucosa is constantly exposed to microbial 
and dietary antigens, all of which are considered to induce 
IgA-secreting plasma cells. Thus, the limited IgA repertoire 
diversity predicted in previous studies (Dunn- Walters et al., 
1997, 2000; Holtmeier et al., 2000; Stoel et al., 2005; Yuvaraj 
et al., 2009) is difficult to reconcile with the broad range of 
intestinal antigens. In this study, we show that in fact each 
individual harbors a private polyclonal and highly diverse IgA 
repertoire. Because we analyzed switched Ig sequences, reper- 
toire analyses performed in this study did not require sorting 
of plasma cells. Instead, RNA was isolated directly from in- 
testinal tissue, and IgA-encoding transcripts were enriched 
during reverse transcription. This approach was equivalent to 
our results obtained when RNA was isolated out of sorted 
plasma cells. The IgA repertoire showed characteristics remi- 
niscent of the T cell receptor repertoire of CD8 + memory 
cells (Naumov et al., 2003), i.e., the IgA repertoire com- 
prised two components which both contributed to the over- 
all IgA repertoire diversity. A first component contained 
clones present at low frequency, whereas a second compo- 
nent encompassed highly expanded clones. Previous studies 
(Dunn-Walters et al., 1997; Holtmeier et al., 2000; Stoel 
et al., 2005; Yuvaraj et al., 2009), which lacked the coverage 
to distinguish these components or wrongly assumed that 
expanded clones are representative of the entire plasma cell 
population, inevitably underestimated the number of B cell 
precursors giving rise to the overall plasma cell population. 
Expanded clones were readily detected in previous studies, 
and clonal expansion has been suggested to occur in the 
periphery (Dunn- Walters et al., 1997, 2000; Holtmeier et al., 
2000) and/or locally in the lamina propria (Husband and 
Gowans, 1978; Yuvaraj et al., 2009). Several observations in 
this study support the former hypothesis. First, expanded 
clones were evenly distributed along the SI but not large 
intestine. Second, SI- and large intestine— specific sorting 
of specificities required chemokine receptor CCR9, and 
third, proliferated plasma cells were randomly distributed 
throughout SI villi. All of these observations suggest that 
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clonal expansion is likely to occur predominantly before 
plasmablast entry into the lamina propria. 

Comparing the IgA sequence repertoire of different mice, 
we observed very little overlap. Antibody heavy chain se- 
quences are generated in a series of recombination events 
using multiple V, D, and J segments and further sequence 
diversification at their junctions (Jung et al., 2006). As a con- 
sequence, the number of possible antibody sequences exceeds 
the number of plasma cells. Still, a comprehensive analysis of 
the B cell repertoire in zebrafish revealed a highly stereo- 
typed usage of VDJ combinations in young fish (Jiang et al., 
2011). Moreover, adult zebrafish also showed signs of con- 
version of their antibody repertoires; i.e., different individu- 
als made the same antibodies (Weinstein et al., 2009; Jiang 
et al., 2011). In fact, we also observed some CDR3 sequences 
shared between different mice, and three of these CDR3 
sequences were highly expanded. At present, we cannot 
reliably calculate the probability to observe such overlap 
in nonconvergent repertoires. Still, most CDR3 sequences 
were unique and not shared between mice. Moreover, the 
similarity of the sequence repertoire was not any different 
comparing mice housed in the same or different cages. This 
observation is reminiscent of the observation made in zebra- 
fish that in older animals, expansion of distinct clones con- 
founded similarities of the original, nonexpanded antibody 
repertoire. Thus, effectively, each mouse harbored an individual 
IgA repertoire. 

In young mice, the IgA repertoire was dominated by 
expanded clones carrying few SMs. We therefore propose 
that early in life, each individual selects and expands a set of 
founder B cells, which best fit the actual intestinal antigen 
load to set up the initial pool of intestinal plasma cells. In- 
triguingly, the clones selected at this early time might also 
serve as basis for the further modulation of the IgA repertoire. 
Frequencies of SMs progressively increased during aging, and 
the IgA repertoire installed after plasma cell depletion was 
largely similar to the repertoire before depletion. Because 
plasma cells are terminally differentiated cells, both observa- 
tions could be explained by the existence of circulating IgA- 
positive B cells. Such cells might be established along with 
the first wave of plasmablasts seeding the SI and considered 
IgA memory B cells. Subsequently, their repeated reentry 
into germinal centers could allow for the progressive accu- 
mulation of SMs and result in the complex clonal populations 
of plasma cell we observed in aged mice. At the same time, 
an IgA memory— like population might enable the rapid recall 
of previously selected specificities, including SMs, after 
plasma cell depletion. Alternatively, the recall of identical 
specificities after plasma cell depletion might be viewed as 
amplification of ongoing immune responses. In fact, intesti- 
nal bacteria are permanently taken up into PPs and carried by 
gut dendritic cells into the draining mLNs (Macpherson and 
Uhr, 2004). Thus, the intestinal immune system is under 
constant stimulation, and GALTs might permanently ac- 
commodate activated B cells of the respective specificities. 
Indeed, recent work has shown that the IgA system lacks 



prime-boost characteristics (Hapfelmeier et al., 2010), an im- 
portant hallmark of memory responses. In any case, our 
observations demonstrate that functionally the IgA system 
can recall previously selected specificities and that an individ- 
ual's private IgA repertoire is fixed in the sense that a limited 
set of expanded clones is maintained and vitally contributes 
to the overall plasma cell pool. Nevertheless, depletion of 
plasma cells has been shown to result in loss of self-reactive Ig 
and to improve the chnical outcome in experimental auto- 
immune models (Neubert et al., 2008), and further repertoire 
studies will be needed to see to what extent plasma cell de- 
pletion can permanently eradicate disease promoting self- 
reactive specificities also in the intestine. 

SHM was triggered by microbiota, which on the one 
hand stimulate the maturation of GALTs and germinal cen- 
ters and on the other hand provide a complex array of anti- 
gens. Both of these functions might contribute to the increase 
of SM frequencies observed after colonization of germ-free 
mice. Moreover, SHM was strictly dependent on T cells and 
the presence of ROR"yt. Besides T cell— dependent germinal 
center reactions, in mice, several alternative pathways of 
plasma cell generation, including peritoneal Bl cells, T cell- 
independent processes, and class switch recombination in the 
lamina propria, have been proposed (Suzuki et al., 2010; 
Cerutti et al., 2011). Our data do not provide strong argu- 
ments to support or dismiss the relevance of such pathways 
in young mice. However, in aged mice, similar to human IgA 
sequences (Barone et al., 2011; Gibbons and Spencer, 2011), 
most V H sequences carried numerous SMs. Thus, in aged 
mice, most IgA-secreting plasma cells seem to be generated 
with the help of T cells. Even though we cannot easily com- 
pare the age of humans and mice, one might speculate that 
the comparison of young mice, which still undergo substan- 
tial modulation of their IgA repertoire, with humans might 
overemphasize the presumed differences between both spe- 
cies. Another prerequisite of SHM seems to be the transcrip- 
tion factor ROR7t, and ROR7t~ / ~ mice did not acquire 
SMs even though these mice possess a polyclonal T cell pool. 
One possible explanation of this observation might empha- 
size the function of GALTs and germinal centers for SHMs. 
Yet we cannot rule out the possibility that the absence of 
lymphoid tissue inducer cells in ROR^t^ mice or other 
defects in these mice might also have more direct effects on 
some pathways of IgA generation. Thus, further experiments 
will be needed to dissect the contribution of germinal cen- 
ter—dependent and —independent pathways of IgA induction 
in young and aged mice as well as in humans. 

In sum, we propose two processes that simultaneously 
shape the mouse intestinal IgA repertoire. On the one hand, 
new clones can be recruited into the intestine. This process 
might reflect the cumulative encounter of IgA-inducing 
stimuli throughout life, which the system cannot appropri- 
ately account for on the basis of expanded clones already rep- 
resented in the repertoire. On the other hand, each individual 
harbors a private set of expanded clones that undergo contin- 
ued T cell— dependent maturation. This process coincides 
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with the establishment of stable intestinal communities, and 
changes in the microbiota might be a major factor driving the 
mutation of expanded clones. However, the host immune 
system also reciprocally affects the composition of the micro- 
biota (Fagarasan et al., 2002; Wei et al., 2011), and such a 
process might rather be regarded as reciprocal shaping of the 
host's IgA repertoire along with its intestinal microbiota. 
Together, these processes account for the dynamic changes 
in the IgA repertoire described herein, offer an explanation 
for the repertoire's unexpected stability after plasma cell 
depletion, and provide a framework to further dissect the 
generation and function of intestinal IgA. 

MATERIALS AND METHODS 

Mice. C57BL/6, CCR9 _/ ~, CD3~ y ~, and ROB^t" 7 " mice (all back- 
crossed to C57BL/6 background for at least eight generations) were bred 
under SPF conditions or germ-free conditions at the central animal facility 
of the Hannover Medical School. Data displayed in Fig. 7 A are based on the 
analysis of ex— germ-free mice previously reported in Pabst et al. (2006). 
PP-free mice were generated by injecting pregnant C57BL/6 mice on ges- 
tational day 14 i.p. with 2 mg lymphotoxin receptor [3— Ig fusion protein 
(Krege et al., 2009). Mice were carefully inspected to lack PPs before inclu- 
sion in the study. To label proliferating cells, mice were injected i.p. with 
3 mg BrdU and thereafter continuously received 80 mg/100 ml BrdU 
with the drinking water. For plasma cell depletion, mice were injected i.p. 
on two consecutive days with 20 j.ig Bortezomib (Velcade). All animal 
experiments were performed in accordance with institutional guide- 
lines and have been approved by the review board of the Hannover 
Medical School and the Niedersaechsische Landesamt fur Verbraucherschutz 
und Lebensmittelsicherheit. 

Mice surgery. Mice were anesthetized using ketamine and xylazine. The 
abdominal cavity was opened, and an SI fragment was exposed. Approxi- 
mately 1x2 mm 2 of intestinal tissue was cut from the an times enteric site, 
and the intestine was closed with suture. The SI was reintroduced into the 
abdominal cavity. For splenectomy, the peritoneum was opened, the spleen 
was located, vessels and ligaments were heat cauterized, and the spleen was 
cut out. The peritoneum was closed with suture, and skin was stapled. 

Cell isolation, flow cytometry, and cell sorting. For isolation of intes- 
tinal lamina propria cells, gut content and PPs were removed, and intestines 
were opened longitudinally and washed with cold PBS. Intestines were 
incubated three times for 15 min in 15 ml HBSS containing 10% FCS and 
2 niM EDTA at 37°C to remove epithelial cells. After each incubation step, 
tubes were shaken vigorously for 10 s, and media containing epithelial cells 
and debris were discarded. The remaining tissue was incubated for 45 min in 
RPMI 1640 with 10% FCS, 0.24 mg/ml collagenase A (Roche), and 40 U/ml 
DNase I (Roche), and tubes were shaken vigorously for 10 s. Cell suspen- 
sions were filtered through a nylon mesh and purified by density gradient 
centrifugation with 40—70% Percoll (GE Healthcare). For flow cytometric 
analysis, cells were stained with IgA-FITC (Invitrogen) and CD138-PE 
(BD). BrdU staining was performed with the APC-BrdU Flow kit (BD) 
according to the manufacturer's instructions. Plasma cells were sorted as 
DAPI"IgA + CD138 + single cells to >95% purity. 

ELISPOT assay. ELISPOT was performed to detect IgA-secreting cells. 
ELISPOT plates (MultiScreen HTS IP 0.45 um; Millipore) were coated 
overnight at 4°C with affinity-purified anti— mouse IgA (5 ug/ml in PBS; 
eBioscience). The next day, plates were washed twice with RPMI 1640 and 
blocked with 10% FCS in RPMI 1640 at 37°C for 2 h. Isolated lamina pro- 
pria cells were added to the wells and incubated overnight at 37°C and 5% 
C0 2 . The next day, cells were discarded, and plates were washed twice with 
deionized H 2 0 and three times with 0.1% Tween in PBS. Detection of IgA 



spots was performed with biotinylated anti-IgA antibody (AbD Serotec) 
followed by peroxidase-labeled streptavidin (Jackson ImmunoResearch 
Laboratories, Inc.) and 3-amino-9-ethylcarbazole (Sigma-Aldrich) as sub- 
strate. Color development was stopped by washing several times with water. 
Plates were scanned, and spots were counted. 

Histology. For preparation of paraffin sections of the SI, a distal piece of the 
SI was washed with PBS and fixed with Bourn's solution (Sigma-Aldrich) for 
5 h at 4°C. After washing three times for 15 min with cold PBS, the SI was 
dehydrated with ascending concentrations of ethanol followed by xylene and 
finally embedded in paraffin. 6-u.m sections were prepared and stained 
according to the manufacturer's instructions with the BrdU Immunohisto- 
chemistry System (EMD) based on Cyanine3-labeled tyramide (PerkinElmer) 
as substrate in combination with IgA staining by an FITC-labeled antibody 
(Invitrogen). Staining was documented using an epifluorescence microscope 
(BX61; Olympus) and AnalysisD software (Soft Imaging System). 

454 sequencing. Three pieces of 5-cm length each were taken from the 
proximal, middle, and distal part of the SI. Gut content and PPs were re- 
moved before pieces were opened longitudinally and washed with cold PBS. 
Intestines were homogenized in TRIZOL reagent (Invitrogen) with an 
Ultra Turrax (all three parts together in 5 ml or each part in 3 ml; IKA). 
From 1 ml of the lysate, RNA was isolated according to the manufacturer's 
instructions. Biopsy samples were lysed in 0.5 ml TRIZOL, and the whole 
lysate was used for RNA isolation. For sequence analysis from sorted cells, 
lamina propria cells were isolated, IgA + CD138 + live (DAPI~) plasma cells 
were sorted, and RNA was isolated with the RNeasy Mini kit (QIAGEN). 
cDNA synthesis was performed with Superscript III (Invitrogen) based 
on a mix of three gene-specific primers for the IgA locus (5'-ATCAG- 
GCAGCCGATTATCAC-3', 5'-TCTCCTTCTGGGCACTCG-3', and 
5'-TGAATGATGCGCCACTGT-3'). To generate template libraries of 
rearranged IgA sequences from IgA + plasma cells, PCR with a primer binding 
in the constant Cot region 5'-CGTATCGCCTCCCTCGCGCCATCAG 
(MID)GAGCTCGTGGGAGTGTCAGTG-3' in combination with a pro- 
miscuous V H primer 5'-CTATGCGCCTTGCCAGCCCGCTCAGGAG- 
GTGCAGCTGCAGGAGTCTGG-3 ' (binding to all V H genes) was 
performed. MID (multiple identifier) specifies a 4-nt sequence used in six 
different variations to identify samples within single lanes. PCR conditions 
were as follows: 95°C, 4 mm; 25 x (94°C, 30 s; 62°C, 30 s; 72°C, 35 s); and 
72°C, 10 min. Amplicons were purified by gel extraction (QIAquick Gel 
Extraction kit; QIAGEN) and quantified by Quant-iT dsDNA HS Assay 
kit (Invitrogen) measured with the Qubit fluorometer (Invitrogen). Ampli- 
cons were prepared with the GS FLX Titanium SV emPCR kit (Lib-A) for 
ultra-deep 454 pyro sequencing on the Genome Sequencer FLX system 
(Roche) as described by the manufacturer. 

Sequence analysis. Sequences were assigned to individual samples ac- 
cording to their MID. Sequences lacking one or both primer sequences or 
shorter than 320 bp were excluded. Sequences were further analyzed 
with ImMunoGeneTics (IMGT) HighV-QUEST (http://www.imgt.org/; 
Brochet et al., 2008; Lefranc et al., 2009), a web portal allowing for the 
analysis of high numbers of sequences. All sequences were compared against 
reference sequences from the IMGT database. Results obtained from IMGT 
were further analyzed with Excel (Microsoft) and VBA {Visual Basic for 
Applications), and only productive sequences were used. Mutation frequencies 
were calculated as the number of mutations divided by the number of all 
nucleotides of the given frame work regions and CDRs. Venn diagrams 
were generated with VennMaster 0.37.4 (2010-04-06). For phylogenetic 
analysis, sequences were aligned with ClustalW 2.0.12 Multiple Sequence 
Alignment tool (Larkin et al., 2007; Goujon et al., 2010). Based on the 
alignment, phylogenetic trees were calculated with ClustalW2 and displayed 
with dendroscope 2.7.4 (Huson et al., 2007). 

Repertoire enumeration and analysis. CDR3 sequences were sorted 
by frequency, and sequences with the same frequency were grouped 
into ranks ranging from single sequences with rank 1 and high abundant 
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sequences given a rank number according to their frequency. To display 
the power law— like distribution, a log-log transformation of a frequency 
versus rank plot was performed. As described in the first section of Results, 
low frequency sequences follow the power law and lay on a straight line, 
whereas high frequency sequences do not. These two components of the 
repertoire can be divided at a critical point, which was calculated as Xc = 
exp{{log N + log a)/— b), where N is the number of different CDR3 
sequences, and a and b are the intercept and slope of the log-log plot. Of 
note, all ranks belonging to the power law component must have a smaller 
rank than Xc. 

Statistical analysis. Statistical analysis was performed with Prism software 
(GraphPad Software). All significant values were determined using one-way 
analysis of variance with Tukey's post-hoc test, for comparing more than 
two groups, or an unpaired two-tailed Student's t test, for comparison of 
two groups. Data are mean + SD. P-values are indicated as follows: *, P < 
0.05; **, P < 0.01; ***, P < 0.001. 
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